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Abstract. The Kovacs protocol, based on the temperature shift experiment originally 
conceived by A.J. Kovacs and applied on glassy polymers jl], is implemented in an exactly 
solvable model with facilitated dynamics. This model is based on interacting fast and slow 
modes represented respectively by spherical spins and harmonic oscillator variables. Due to 
this fundamental property and to slow dynamics, the model reproduces the characteristic non- 
monotonic evolution known as the "Kovacs effect", observed in polymers, spin glasses, in 
granular materials and models of molecular liquids, when similar experimental protocols are 
implemented. 



1. Introduction 

Glassy systems, being in an out-equilibrium condition, have properties which depend on their 
history. This is the 'memory' of glasses. This property can manifest itself in striking ways, when 
specially devised experiments are made. One example is given by negative temperature cycles in 
spin glasses where the ac susceptibility, depending on both frequency and the age of the system, 
recovers the exact value it had before the negative temperature jump. A memory effect which 
shows up in a one-time observable, when a specific experimental protocol is implemented, is the 
so called "Kovacs effect" pQ. This effect has been the subject of a variety of recent investigations 
(21 EH IU El El El- ^he characteristic non-monotonic evolution of the observable under examination 
(the volume in the original Kovacs' experiment), with the other thermodynamic variables held 
constant, shows clearly that a non-equilibrium state of the system cannot be fully characterized 
only by the (time-dependent) values of thermodynamic variables, but that further inner variables 
are needed to give a full description of the non-equilibrium state of the system. The memory in 
this case consists in these internal variables keeping track of the history of the system. 

The purpose of this paper is to use a specific model for fragile glass to implement the protocol 
and get some insight into the Kovacs effect. We show that in spite of its simplicity, this 
model captures the phenomenology of the Kovacs effect and allows in specific regimes to obtain 
analytical expressions for the evolution of the variable of interest. This paper is organized as 
follows: in Section II we review the experimental protocol generating the effect, in Sections III 
and IV we introduce our model and use it to implement the protocol, in section V we draw out 
of this model some analytical results with final conclusions. An appendix collects all terms and 
coefficients employed in the main text. 



2. Kovacs protocol 

The experimental protocol, as originally designed by A. J. Kovacs in the '60s fQ, consists of 
three main steps: 

1 st step The system is equilibrated at a given high temperature Tj. 

2 nd step At time t = the system is quenched to a lower temperature 7], close to or below 
the glass transition temperature, and it is let to evolve a period t a . One then follows the 
evolution of the the proper thermodynamic variable (in the original Kovacs experiment 
this was the volume V(t) of a sample of polyvinyl acetate, in our model it will be the 
"magnetization" M\(t)). 

3 rd step After the time t a , the volume, or other corresponding observable, has reached a 
value equal, by definition of t a , to the equilibrium value corresponding to an intermediate 
temperature Tf (T\ < Tf < Tj), i.e. such that V^ita) = V^ q . At this time, the bath 
temperature is switched to Tf. 

The pressure (or corresponding variable) is kept constant throughout the whole experiment. 

Naively one would expect the observable under consideration, after the third step, to remain 
constant since it already has (at time t = ) its equilibrium value. But the system has not 
equilibrated yet and so the observable goes through a non monotonic evolution before relaxing 
back to its equilibrium value, showing a characteristic hump whose maximum increases with the 
magnitude of the final jump of temperature Tf — Ti and occurs at a time which decreases with 
increasing Tf — T\. 

We want to implement this protocol on a model for both strong and fragile glass first 
introduced in the Harmonic Oscillators-Spherical Spins model (HOSS). This model is based 
on interacting fast and slow modes, this property turns out to be necessary for the memory 
effect, object of this paper, to occur. 



3. The Harmonic Oscillator- Spherical Spin Model 

The HOSS model contains a set of N spins Si locally coupled to a set of N harmonic oscillator 
Xi according to the following hamiltonian: 

- K 

n = E(y x * 2 - Hxi ~ JXiSi ~ L5 *) W 
i=i 

The spins have no fixed length but satisfy the spherical constraint: YliLi = N. The spin 
variables are assumed to relax on a much shorter time scale than the harmonic oscillator 
variables, so the oscillator variables are the slow modes and on their dynamical evolution the 
fast spin modes act just as noise. One can then integrate out the spin variables to obtain 
the following effective Hamiltonian for the oscillators (for details see jS], explicit expressions of 
undefined terms appearing in all equations hereafter are reported in the Appendix): 
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which depends on the temperature and on the first and second moment of the oscillator variables, 
namely: 

N N 
i=l i=l 



These variables encode the dynamics of the system which is implemented through a Monte Carlo 
parallel update of the oscillator variables: 



Xi -> Xi + n/VN (4) 

The variables r{ are normally distributed with zero mean value and variance a 2 . The update is 
accepted according to the Metropolis acceptance rule applied to the variation 5e of the energy 
of the oscillator variables, which is determined by 7~i e ff and, in the limit of large N, is given by: 

| = K T {Mi,M 2 ) 5M2 _ ffT(Mi)Mam _ (5) 

This simple model turns out to have a slow dynamics and can be solved analytically. 
Following [S] one can derive the dynamical equations for M\ and M 2 



Mi 



H T (M 1 ,M 2 ) 



K T (M 1 ,M 2 ) 

Mo = 

K T (M 1 ,M 2 ) 



f T (M l ,M 2 ) (6) 
J T (Mi, M 2 ) + Ht{M\, M 2 )M\ 



The stationary solutions of these equations coincide with the saddle point of the partition 
function of the whole system at equilibrium at temperature T and are given by: 

= H T (M U M 2 ) Ht 

K t (Mi,M 2 ) K T v ; 

T T 

1 K T (M 1 ,M 2 ) K T 

with barred variables from now on indicating their equilibrium values. 
3.1. Strong and Fragile Glasses with the HOSS model 

In spite of its simplicity, the HOSS model allows to describe both strong and fragile glasses, 
characterized respectively by an Arrhenius or a Vogel-Fulcher law in the relaxation time. The 
following constraint on the configurations space is applied: 



1 1 2 



M 2 - Ml - M > (8) 



When Mq = there exists a single global minimum in the configurations space of the oscillators, 
therefore the role of the constraint with Mq > is to avoid the existence of a "crystalline state" 
and to introduce a finite transition temperature. The stationary solutions for the dynamics with 
this constraint are given by: 

= H T (M U M 2 ) _H T 
K T {M U M 2 ) K T 

T<T k 

The temperature is determined by the further condition: 



M 2 - Ml 




T k = M K Tk (M^ k , M 2 Tfc ) = M K Tk . 



(10) 



This is the highest temperature at which the constraint is fulfilled, for smaller temperatures the 
system relaxes to equilibrium configurations which fulfill the constraint. For T > therefore 
the dynamics is not affected by the constraint. For T < Tk the system eventually reaches 
a configuration which fulfills the constraint, when this happens it gets trapped for ever in 
such a configuration. This is equivalent to having a "Kauzmann-like" transition, occurring at 
T = Tk with vanishing configuration entropy, meaning the system gets stuck forever in one single 
configuration fulfilling the constraint (see also: |5|). 

When there is no constraint, i.e. when Mq = 0, then Tk = 0, if the Monte Carlo updates 
are done with Gaussian variables with constant variance a 2 , this model is characterized by an 
Arrhenius relaxation law: 

r eq ~ e~ (11) 

in so resembling the relaxation properties of strong glasses. 

The HOSS model with constraint strictly positive (Mo > 0) can easily be extended to describe 
fragile glasses by further introducing in the variance of the Monte Carlo update, the following 
dependence on the dynamics: 

a 2 = 8(M 2 - M 1 2 )(M 2 - Ml - M ) -7 (12) 

In this case the relaxation time turns out to follow the generalized Vogel-Fulcher law: 

Teq ~ e ^/(T-nr (13) 

The parameter 7 is introduced to make the best Vogel-Fulcher type fit for the relaxation time in 
experiments, making this model valid for a wide range of fragile glasses. When the temperature 
approaches the value Tk defined by ()10j) . from above, the system relaxes towards configurations 
close to the ones fulfilling the constraint. The variance a 2 then tends to diverge, the updates 
become large and so unfavorable, meaning that almost every update of the oscillator variables 
is refused. This produces the diverging relaxation time following the Vogel-Fulcher law of Eq. 

(USD. 



4. Kovacs effect in the HOSS model 

We implement the Kovacs protocol in the model above introduced for a fragile glass. The system 
is prepared at a temperature T L and quenched to a region of temperature close to the Ty., i.e. 
T[ > Tfc. Solving numerically Eqs. (JHJ) we determine the evolution of the system in both step 2 

and 3 of the protocol. In step 2 the time t a , at which M 1 1 (t a ) = M x 1 , is calculated so that: 

M?(fi) = M? (14) 
M?(t+) = M?{t a ) 

The evolution of the fractional " magnetization" : 

^um^MiB^t (15) 

after step 3 (t > t a ) for different values of 7] is reported in Figs. ^ and [21 respectively for 
7 = 1 and 7 = 2. The magnetic field H is kept constant at the value H = 0.1. In all the 
implementations of the protocol we use the values J=K=1, L=0.1 and Mo=5 for the other 
parameters of the model. This choice for the parameters and the value H=0.1 for the magnetic 
field fix (through Eqs. @ and 1)101) ) the Kauzmann temperature at the value T^=4. 00248. 




Figure 1. Fragile glass with 7=1 The 
Kovacs protocol is implemented with a 
quench of the system from temperature 
Tj=10 to Ti, and final switch (at t=t£) to 
the intermediate temperature Ty=4.3. The 
continuous lines, starting from the lowest, 
refer to 7}=4. 005, 4.05, 4.15, the dashed line 
refers to condition T\ = Tf, ( simple aging 
with no final temperature shift). 



Figure 2. Fragile glass with 7=2. 

The Kovacs protocol is implemented with 
a quench of the system from temperature 
Tj=10 to Ti, and final switch (at t=t£) to 
the intermediate temperature Tj=4.3. The 
continuous lines, starting from the lowest, 
refer to T,=4.005, 4.05, 4.15, 4.25, the dashed 
line refers to condition T\=T$ (simple aging 
with no final temperature shift). 



Since the equilibrium value of M\ decreases with increasing temperature (as opposed to what 
happens for instance with the volume) we observe a reversed 'Kovacs hump'. The curves keep 
the same properties typical of the Kovacs effect, the minima occur at a time which decreases 
and have a depth that increases with increasing magnitude of the final switch of temperature. 
As expected, since increasing 7 corresponds to further slowing the dynamics, the effect shows 
on a longer time scale in the case 7 = 2 as compared to 7 = 1. 

rp 

Actually, since in the last step of the protocol: M\(t = t a ) = f and fa J Mi, M 2 ) is always 
positive, from the first of Eqs. © one soon realizes that the hump for this model can be either 
positive or negative, depending on the sign of the term: 

H Tf (M?,M 2 ) mT} 
K Tf (M^,M 2 ) 1 

at t = t^. This term is zero when Mi = M^ s ,M 2 = , so one would expect M 2 (t = 

_ rp 

t£) = M 2 to be the border value determining the positivity or negativity of the hump. Since 

{M 1 ,M2) decreases with increasing M 2 while K Tf (M 1 ,M 2 ) increases, it follows that the 
condition for a positive hump is: 

M 2 (t = tt) < M 2 f (17) 

For shifts of temperature in a wide range close to the transition temperature T/%, where the 
dynamics is slower and the effect is expected to show up significantly on a long time scale, the 

_ rp 

condition M 2 (t = t a ) > M 2 f is always fulfilled and therefore a negative hump is expected. 



(16) 



5. Analytical solution in the long-time regime 

In the previous Section we have shown, through a numerical solution of the dynamics, that the 
HOSS model reproduces the phenomenology of the Kovacs effect, showing the same qualitative 



properties of the Kovacs hump as obtained in experiments (see for ex. PQ2j), in some other 
models with facilitated or kinetically constrained dynamics E] and in other different models 
PI IS 13 IE] • I n this section we show that, by carefully choosing the working conditions in which 
the protocol is implemented, our model provides with an analytical solution for the evolution of 
the variable of interest. 



5.1. Auxiliary variables 

In order to ease calculations, as done in 015] it is convenient to introduce the following variables: 



Ml 



H T (M U M 2 ) 
K T (M 1 ,M 2 ) 



Mi 



(18) 



H2 = M 2 - M{ 



Mn 



for which the dynamical equations read: 



Ml 
h 



JQ T (M 1 ,M 2 )I T (M 1 ,M 2 ) - (l + DQ T (M 1 ,M 2 ))fJLifT(M 1 ,M 2 ) 
+ 2f4f T (M 1 ,M 2 ) 



(19) 



2I T (M 1 ,M 2 ) 2 



K T {M U M 2 ) 



We will choose to implement steps 2 and 3 of the protocol in a range of temperature very 
close to the Kauzmann temperature T^. As exhaustively shown in [HJ^J in the long time regime 
the variable fj> 2 (t) decays logarithmically to its equilibrium value which is small for T ~ Tk- So, 
if t a is very large, the value of the variable ^ 2 (t), which is continuous at the jump, will be small 
enough to fulfill the condition for which the following equation is shown to be valid [Hj: 



dm 



d(5fi 2 ) 



{l + Q T {M 1 ,M 2 )D) 



(p2 + 6fl 2 ) 
#M2 



"Mi 



JQ T {M U M 2 )T 
2 (M + fl 2 ) 



(20) 



where now the variable Sfx 2 (t) = fJ- 2 (t) — p, 2 is used and barred variables always refer to 
equilibrium condition. Of course choosing 7} close to and waiting a long time t a so that 
the system approaches equilibrium, allows only small temperature shifts for the final step of the 
protocol, meaning that also Tf will be close to TV All the coefficients which appear in equation 
(|2U|) (see Appendix for complete expressions) in the regime chosen, can be assumed constant 
and equal to their equilibrium values with a very good approximation. The equation can then 
be easily integrated to give: 



Mi(^M2 



exp 



2^1(7,7,7 + 1, 



J£2_' 
5fl2 • 



7(5/22 p/A Q 



1-4 b q - Cq dzexp 

J 5^ 



2^(7,7,7 + l,-f) 
1^/A Q 



where the superscript + indicates t = t+ and 2-^1 the hypergeometric function. This expression 
simplifies in cases 7=1, 3/2, 2. All these solutions and relative coefficients are reported in the 
appendix, here we limit ourselves to the case 7 = 1 which corresponds to ordinary Vogel-Fulcher 
relaxation law. In this case the solution is: 



MiW 



*M2(*) 



where: 



5/i 2 (t) + fx 2 



dz 



Mi^ 



6f4 



M2 



5/4 



fa 

M2 



Cq 



5fM 2 (t) 



Stq 



z + rj 



x a+1 2*1(1 



a, —a, 2 



a, 



dz 



"-) 



z + M2 



.fa 

M2 



(21) 
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Figure 3. Comparison between numerical solution (continuous lines) for the Kovacs' curves and 
the approximate analytical solution at short-intermediate (dot-dashed line) and intermediate- 
long time (dashed line). The protocol is implemented between = 10 and Tf = 4.018. The 
curves starting from the lowest refer to T t = 4.005, 4.008. (H = 0.1, T k = 4.00248) 

One can then expand the variable of interest M\(t) in terms of jX\ and 5fi2 and obtain the 
following expression for the Kovacs curves: 

5M X (t) = A\ f (mJ> , M 2 T/ ) ( Ml (t) - M +) + A\ f (M? f , ) (fy 2 (t) - «5 M +) (22) 

where the coefficients are approximately constant in the regime chosen and can be evaluated at 
equilibrium. 

5.2. Short and intermediate t — t a 

For small t — t a , a linear approximation for the variable 5fi2, with slope given by the second 
equation of the set (|19|) evaluated at t = t+, turns out to be very good. Inserting this expression 
in Eq. 1)21(1 to get ^\{t) and then in Eq.(|22|) a good approximation of the first part of the hump 
for small and intermediate t — t a is obtained, as shown in Fig. O 

5. 3. Intermediate and long t — t a 

When t — t a is very large, we can use Eq. (|21j) and the pre-asymptotic approximation for ^{t) 
(see: 0) 

fi 2 (t) = (log{t/t ) + \ log(log(t/to))) (23) 

Inserting this expression in Eq. (j2*Tj) to get Hi(t) and then in Eq. (|2*2*|) . a good approximation 
for the hump and the tail of the Kovacs curves is obtained. In Fig. El we show the agreement 
between the analytical expression so obtained and the numerical solution. 

6. Conclusion 

We have shown that a simple mode with constrained dynamics like the HOSS model, is rich 
enough to reproduce the Kovacs memory effect, even allowing to obtain analytical expression 
for the Kovacs hump in a long time regime. The Kovacs effect is observed in many experiments 
and models, showing common qualitative properties which we have found to be shared also by 
the model analyzed in this paper. The quantitative properties depend on the particular system 
or model analyzed. 



As far as it concerns the HOSS model, it turns out that for the slow modes, i.e. the 
oscillator variables, fixing the overall average value, the magnetization Mi, does not prevent 
the existence of memory encoded in the variable M2, which keeps track of the history of the 
system. The equilibrium value of M 2 increases with temperature while the equilibrium value 
of Mi decreases with increasing temperature. Therefore after the final switch of temperature, 
since M 2 (t a ) > M 2 , the variable M2 has a value corresponding to an equilibrium condition at a 
higher temperature (memory of the initial state at temperature Tj) so driving the system towards 
a condition corresponding to a higher temperature, i.e. smaller values of Mi, determining the 
hump. 

It is important to stress that a fundamental ingredient in the HOSS model, besides the slow 
dynamics which originates from the Monte Carlo parallel update, is the interaction between slow 
and fast modes. Due to this interaction the equilibrium configurations of the oscillator variables 
at a given temperature are determined by both M2 and Mi, the first and second moment of their 
distribution, whose dynamical evolution is interdependent. When such interaction is turned off 
(by setting J = ) essentially only one variable is sufficient to describe the equilibrium con- 
figurations and the dynamics of the system, and the memory effect is lost. In this respect this 
model constitutes an improvement to the so-called oscillator model ^3] within which such mem- 
ory effect cannot be reproduced. In the present model one can also study temperature cycle 
experiments of the type carried out in spin glasses (see. ^U), leaving room for further research. 
More details on this subject can be found in Ref. [To] . 

G.A. and L.L. gratefully acknowledge the European network DYGLAGEMEM for financial 
support. 



Appendix A. 

In this Appendix we report all the explicit expressions for terms appearing in the text. In Eqs. 
©, © and © we have: 



w T (M 1 ,M 2 ) = ^/J 2 M 2 + 2JLM 1 + L 2 + T 2 /4 

J 2 



K T (M 1 ,M 2 ) = K 



w T {M 1 ,M 2 )+T/2 
H T (Mi,M 2 ) = H + — x - 

/t(Mi,M 2 ) = - ^^l^- Erfc[a T {M u M 2 )] • exp [4(M 1; M 2 ) - 4(Mi,M 2 )] 
7 T (Mi,M 2 ) = ° Kt ^ X ' M2) Erfc[a T {M 1 ,M 2 )\ + ( ^ - K t (M 1: M 2 )wt(M u M 2 ) ) / T (Mi,M 2 ) 



2 



where: 



a T (Mi,M 2 ) = 



y Sw T {M l ,M 2 ) 

d T (M 1 ,M 2 ) _ 2K T (M 1 ,M 2 )w T (M 1 ,M 2 ) 
a T (M 1 ,M 2 ) ~ T 



In Eqs. CETJ), JTJJ), ((201), (EU, (ED) and (|22j): 



A, 



a Kt 



Q t (M 1 ,M 2 ) 
P T (M 1 ,M 2 ) 



D = JH + LK = JH T + LK T 
J 2 D 



2 



K^w T (w T + T/2) 2 

j\M 2 -M 2 ) 
2K t wt(w t + T/2) 2 
K Tk {K-K Tk ){l + DQ Tk +P Tk ) 
(K-K Tk ){l + DQ Tk )-K Tk Q Tk 

Fi{a,b,c,z) = — > —— r r 

r ( a ) r W^ r(c + n) n\ 

Aq = l + QT f (Mp,M^ f )D = l + Q Tf D 



Cq 



2^1(7, 7,7 + 1,-^) 
JQ Tf (M? , M 2 T/ )T f JQ Ts T f 



Bl = e W [A Q ^ 



2(M + /2 2 ) 2(M + jS 2 ) 



M \ (w T + T/2)K T 

A T (M U M 2 ) - Mi{JMi+l + [wt+t/2)Kt) 

A 2 T {M l ,M 2 ) = 2MiA 1 T {M l ,M 2 ) 

1 + DQ T 



to 



87 1 + DQ T + P T 
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